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I. INTRODUCTION 



After the first experimental realization of Bose- 
Einstein condensation in dilute atomic gases [l| the field 
of ultracold gases has become a rapidly gxowingfield of 
research (for reviews see, for example, @, H, 0, [SIS 13 • I n 
the last years a considerable amount of experimental and 
theoretical work has focused on Fermi gases. With re- 
spect to Bose gases Fermi systems exhibit important dif- 
ferences which are the consequence of quantum statistics 
and of the role of interactions. A first important differ- 
ence is that, at low temperature, dilute Fermi gases occu- 
pying a single spin state practically do not interact since 
s-wave scattering is suppressed by the Pauli principle. 
This provides a unique opportunity for an almost perfect 
realization of the ideal Fermi gas, with useful applications 
to precision measurements and quantum information pro- 
cesses. A second important difference is that superfluid- 
ity in Fermi gases is the result of non trivial many-body 
mechanisms where interactions play a crucial role giving 
rise to pairing effects. The resulting many-body state 
exhibits a rich variety of features, depending on the sign 
and the value of the scattering length characterizing the 
interaction between atoms belonging to different atomic 
species. A particularly interesting configuration is the 
so called unitary limit where the scattering length takes 
an infinite value. At unitarity the system is found to be 
particularly stable and to exhibit clean manifestations of 
superfluidity. The possibility of tuning the value of the 
scattering length profiting of the existence of Feshbach 
resonances and the rich variety of trapping configura- 
tions, both of magnetic and optical nature, are making 
the study of ultracold Fermi gases a rich subject of re- 
search with many stimulating opportunities from both 
the experimental and theoretical side. 

The purpose of this paper is to review some of the 
dynamic and superfluid features exhibited by ultracold 
Fermi gases with special emphasis on the effects of the 
external confinement which will be assumed in most cases 
of harmonic shape. After introducing the main features 
of the ideal Fermi gas in a harmonic trap (Sect. 2) we 
will discuss the role of interactions and the general be- 
havior exhibited by an interacting Fermi gas along the 
BEC-BCS crossover (Sect. 3). We will then apply the 
many body results of Sect. 3 to trapped configurations 
(Sect. 4). Sects. 5 and 6 will be devoted to the study of 
superfluid effects, both concerning the dynamic behavior 
(expansion and collective oscillations) and the rotational 
properties (Sect. 6). 



II. IDEAL FERMI GAS IN HARMONIC TRAP 

The ideal Fermi gas represents a natural starting point 
for discussing the physics of dilute Fermi gases. In some 
cases the role of interactions can be in fact neglected (this 
happens for spin polarized gases where the interaction is 
strongly suppressed by the antisymmetrization require- 
ment). Viceversa when the role of interactions becomes 
crucial the comparison with the predictions of the ideal 
gas points out explicitly new interesting features exhib- 
ited by Fermi gases. 

The ideal Fermi gas in the harmonic potential 
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is a subject with many applications in different fields of 
physics, ranging from nuclear physics to quantum dots 
and is widely discussed in the lterature. For this reason 
we will only focus on the most relevant features of the 
model, emphasizing the large N behavior where the mo- 
tion of the gas can be described in semiclassical terms. 
In the large N limit many single particle states are in 
fact occupied and the role of the Heisenberg uncertainty 
principle can be safely ignored for most properites of the 
system. The simplest way to introduce the semiclassical 
description is to use a local density approximation for the 
Fermi distribution function of a given spin species: 
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where /i is the chemical potential fixed by the normaliza- 
tion condition 
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and N a is the number of atoms of the given spin species. 
Equation ^ explicitly shows that, although incompat- 
ible with the Heisenberg indetermination relation, the 
semiclassical approach accounts for the Pauli exclusion 
principle which represents the most relevant feature ex- 
hibited by Fermi gases. In Eq. ([3]) we have introduced 
the single particle density of states 
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where e sp (r,p) = p 2 /2m + V/ lo (r) is the classical energy. 
The density of states depends on the dimensionality as 
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well as on the actual form of the confining potential, 
the harmonic trapping potential |T]) we find 
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where uJho = (^x^y^z) 1 ^ 3 is the geometrical average 
of the three trapping frequencies. The resulting en- 
ergy dependence differs from the one of the uniform 3D 
gas where the density of states takes the value g(E) = 
\/E~V(m 3 / 2 /(y/2h 3 Tr 2 )). The difference has its physical 
origin in the suppression of states in phase space due to 
the spatial confinement produced by the oscillator poten- 
tial. 

In terms of the density of states one can easily calculate 
the relevant thermodynamic functions. For example the 
energy of the gas is given by the expression 
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energy. At T = one has E(0) = (3/A)N a E F . In Fig.l 
we show how the release energy ^ varies as a function of 
the temperature. The comparison with the prediction of 
the classical gas (dashed line) and of the ideal Bose gas 
reveals explicitly the effects of quantum statistics on this 
measurable quantity. These effects were clearly demon- 
strated in the experiment of [H, Q which are reported in 
Fig.2. 

The Fermi energy ([7]) can be used to define typical 
length and momentum scales characterizing the Fermi 
distribution in coordinate and momentum space respec- 
tively: 
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where R x , R y and R z are the widths of the density 
distribution at zero temperature which can be directly 
calculated integrating the T = distribution function 
/(r, p) = 9(e(r, p) — Ep) in momentum space: 



At T — the chemical potential defines the Fermi energy 
(fi(T = 0) = Ep) and the normalization condition yields 
the result 



Ep = k B T F = {QN^huJho 



(T) 



which fixes an important energy (and tempertaure) 
scale in the problem, analog to expression Ep = 
(H 2 /2m)(6ir 2 n (T ) 2 / 3 holding for uniform matter where n a 
is the density of a single spin component. 

It is worth noticing that the Fermi energy ([7]) has the 
same dependence on the number of atoms and on the 
oscillator frequency Who as the critical temperature for 
Bose-Einstein condensation, given by the well know for- 
mula kgT c = 0.94/iaj/joiV 1 / 3 . This is not a surprise since 
in a gas the effects of quantum degeneracy become im- 
portant when the thermal wave length (2ttH 2 / (inksT) 1 ' 2 
is of the order of the average distance n _1 / 3 (0) between 
atoms where n(0) is the density in the center of the 
trap. Using a classical Gaussian distribution to pro- 
vide an estimate of the density of the gas one finds 
n-V 3 (0) - iV- 1 / 3 (fc B T/TOw^) 1 / 2 , so that the scale of 
temperatures relevant for observing quantum effects is 
given by fkUhoN 1 / 3 /Hb in both Fermi and Bose gases. 
In typical experiments the value of Tp corresponds to 
microkelvin or fractions of microkelvins. It is however 
worth noticing that, differently from Tbec, the Fermi 
energy does not define the critical temperature of a phase 
transition, but just a crossover characterizing the onset 
of quantum degeneracy phenomena. The occurrence of a 
phase transition in a Fermi gas can be only the result of 
interaction effects. 

An important quantity to calculate is also the release 
energy E re i given by the energy of the gas after switch- 
ing off the confining potential. A consequence of the 
equipartition theorem applied to the ideal gas trapped 
by the harmonic potential is that the release energy is 
always equal to E Te i — E(T)/2 where E(T) is the total 
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while the Fermi momentum pp fixes the width of the 
corresponding momentum distribution 
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obtained by integrating the T = distribution func- 
tion in coordinate space. Equations ((9]) and (fT0|) hold 
for positive values of their arguments and are often re- 
ferred to as Thomas-Fermi distributions. Equation (JTUJ) 
is the analogue of the most familiar momentum distribu- 
tion 3iV (T /(47rp'p)9(l— p 2 jpp) characterizing the uniform 
Fermi gas. The smoothing of n (p) with respect to the 
uniform case is the consequence of the harmonic trap- 
ping. Notice that the value of pp defined above coincides 
with the Fermi momentum 
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of a uniform gas evaluated in the center of trap. Using 
Eqs.([7]) and © one easiliy finds the useful expressions 



and 
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for the Thomas-Fermi radii and for the Fermi momen- 
tum, respectively, where a^ = \J h/mtuho is the average 
oscillator length. It is worth comparing equations (l9l fT2|) 
with the analogous results holding for a trapped Bose- 
Einstein condensed gas in the Thomas- Fermi limit 0]. 
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For example the density distribution of the BEC config- 
uration is given by the expression 



(cr =1). We will consider the grand canonical many-body 
Hamiltonian 
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with the Thomas-Fermi radii given by Ri — 
aho{^Na/ ahoY^^ho/^i where a is the s-wave scatter- 
ing length characterizing the interaction between bosons. 
The density profiles © and (fl"4"l) do not look dramati- 
cally different. In both cases the radius of the atomic 
cloud increases with the number of atoms although the 
explicit dependence is slightly different (N 1 / 5 for bosons 
and iV 1 / 6 for fermions). Notice, however, that the form of 
the density profiles has a deeply different physical origin 
in the two cases. For bosons it is fixed by the two-body 
interaction, while in the Fermi case it is determined by 
quantum pressure effects. 

In momentum space the Bose and Fermi distributions 
instead differ in a profound way. First, as a consequence 
of the semiclassical picture, the momentum distribution 
of the Fermi gas is isotropic even if the trapping poten- 
tial is deformed, differently from what happens in the 
BEC case. Second, the momentum width of a trapped 
Bose-Einstein condensed gas scales like 1/R and hence 
decreases by increasing N while, according to eqs.([B}- 
(fl"3f the momentum width of a trapped Fermi gas scales 
like R and hence increases with the number of atoms. 
The different behavior reflects the fact that the Heisen- 
berg uncertainty inequality is close to an identity for a 
Bose-Einstein condensate, while in a Fermi gas one has 
APAR ~ Ep/ivho > h. 

Despite its simplicity the ideal Fermi gas can exhibit 
non trivial features. This is the case if one consider 
the combined presence of harmonic plus periodic poten- 
tials which gives rise to Bloch oscillations [l(| and to 
insulating phenomena [T7J ] of relatively easy experimen- 
tal access. Furthermore the effects of Fermi statistics of 
the ideal Fermi gas show up in a non trivial way in an 
anti-bunching behaviour exhibited by the pair correla- 
tion function which has been the object of recent mea- 
surements Ha. Hal. 



III. ROLE OF INTERACTIONS: THE BEC-BCS 
CROSSOVER 

The ideal gas model presented in Sect. 2 provides a 
good description of a cold spin polarized Fermi gas. In 
this case interactions are in fact strongly inhibited by the 
Pauli exclusion principle. When atoms occupy different 
spin states interactions instead deeply affect the solution 
of the many-body problem. This is particularly true at 
very low temperature where they give rise to pairing ef- 
fects bringing the system into the superfluid phase. 

Let us consider the simplest case of a two-component 
system occupying two different spin states hereafter 
called, for simplicity, spin-up (a =f) and spin-down 



h = E / dr *^ r ) (-^r + Vext{r) fI ) * ff(r) 

+ ydrdr'y(r-r')#j(r)*{(r')*|(r')* T ( r ) > ( 15 ) 

written in second quantization where \x is the chemical 
potential and the field operators obey the fermionic anti- 
commutation relations {^ CT (r), & a ,(r')} = 5 ay(J /5(r — r'). 
The external potential V ext and the two-body potential 
V account, respectively, for the external confinement and 
for the interaction between atoms of different spin. The 
number of atoms, fixed by the normalization condition 
J dr < \^/ a (r)\ 2 >= N a , can in general be different for 
the two spin species. In this section we will consider the 
unpolarized case Nf = N± = N/2 in the absence of ex- 
ternal trapping (V ex t = 0). The inclusion of harmonic 
trapping will be discussed in the next section. 

We are interested in dilute gases where the range of 
the interatomic potential is much smaller than the in- 
terparticle distance. Furthermore, we assume that the 
temperature of the system is sufficiently small so that 
only collisions in the s-wave channel are important. Un- 
der these conditions one expects that interaction effects 
will be governed by a single parameter: the s-wave scat- 
tering length a. In this regard one should recall that the 
gaseous phase corresponds to a metastable solution of the 
many-body problem, the true ground state being in gen- 
eral a crystal configuration where the microscopic details 
of correlations are important (to simplify the notation the 
metastable solution will be often also called the ground 
state). The above considerations explain why in Eq. fPo]) 
we have ignored the interaction between atoms occupy- 
ing the same spin state which, in most cases, is expected 
to give rise only to minor corrections, due to the quench- 
ing effect produced by the Pauli principle. Of course this 
picture can change considerably in the presence of p- wave 
resonances. 

In order to investigate in a more tractable way the 
effects of quantum correlations and better understand 
the role played by the scattering length, it is convenient 
to replace the microscopic potential V with an effective 
potential V e t / . Different types of effective potentials can 
be considered. In many applications one introduces the 
regularized zero-range pseudo-potential [l4| . 
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(16) 



where the coupling constant g is related to the s-wave 
scattering length a characterizing the interaction be- 
tween two atoms of different spin by the relationship 
g = ATrh 2 a/m. The regularization accounted for by the 
term (d/dr)r permits to cure the ultraviolet divergen- 
cies in the solution of the many-body Schrodinger equa- 
tion that arise from the vanishing range of the pseudo- 
potential. In general, the regularization is crucial to solve 
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the many-body problem beyond lowest order perturba- 
tion theory, as happens, for example, in the BCS theory 
of superfluidity. 

In this work we will often discuss the predictions of 
quantum Monte Carlo simulations. In these numeri- 
cal approaches the use of the effective potential (fl6|) is 
hard to implement and one must resort to a different 
effective interatomic potential. A convenient choice is 
the attractive square- well potential defined as follows: 
V e ff(r) = —Vq for r < R and V e ff(r) = other- 
wise [other choices for V e ff have also been considered 
in the literature (l5j . The interaction range Ro must be 
taken much smaller than the inverse Fermi wave vector, 
kpRo *C 1, in order to ensure that all the many-body 
properties of the system will not depend on its value. 
The depth Vq of the potential is instead varied so as to 
reproduce the actual value of the scattering length ac- 
cording to the relation a = Rq[1 — tan(K R ) / (K R )], 
where K$ = mVo/ti 2 . Notice that a becomes infinite for 
every new two-body bound state entering the well. 

The above approaches permit to describe the many- 
body features uniquely in terms of the scattering length 
a. These schemes are adequate if the scatterng length is 
the only relevant interaction parameter and the two-body 
scattering amplitude, given by the expansion 



a- 1 + ik + R*k 2 y ' 

can be safely evaluated by keeping only the first two 
terms in the denominator and neglecting the term con- 
taining the effective range R* of the interatomic poten- 
tial. When R* becomes of the order of the inverse Fermi 
wave vector, as happens in the case of narrow Feshbach 
resonances, more complex effective potentials should be 
introduced in the solution of the many-body problem 

In this section we will focus the discussion on the case 
of uniform systems (V ex t = 0) where exact solutions of 
the many body problem are available in some important 
limiting cases. A first example is the dilute repulsive 
gas. In this case interactions can be treated through the 
effective potential p^|) with a positive scattering length 
a. Standard perturbation theory can be applied with 
the small parameter kp a <C 1 expressing the diluteness 
condition for the gas. The expansion of the energy per 
particle up to terms quadratic in the dimensionless pa- 
rameter kpa then yields, at T = 0, the result [ItJ 

E 3 ^ (, 10 , 4(ll-21og2) ,, \ 
N = l EF \! + 9^ a+ 21** {kFa) -) 

(18) 

where Ep is the Fermi energy. The above result is uni- 
versal as it holds for any interatomic potential with a 
sufficiently small effective range R* such that nR* 3 <C I. 
Higher order terms in the expansion (jT5J) will not de- 
pend only on the scattering length a, but will contain an 
explicit dependence on the details of the potential. In 
the case of purely repulsive potentials, such as the hard- 
sphere model, the expansion in Eq. (| 1 8|) corresponds to 



the energy of the "true" ground state of the system. For 
more realistic potentials with an attractive tail, the above 
result describes instead the metastable gas-like state of 
repulsive atoms. This distinction is particularly impor- 
tant in the presence of bound states in the two-body 
problem, as more stable many-body configurations could 
consist of a gas of dimers (see discussion below) with the 
same positive value of a. It is finally worth noticing that 
the repulsive gas described by the equation of state (fT8|) 
does not exhibit superfluidity. 

A second important case that can be solved exactly 
is the dilute Fermi gas interacting with negative scatter- 
ing length (feflal *C 1). In this limit, hereafter called 
the BCS limit, the many-body problem can be solved 
both at T — and at finite temperature and corresponds 
to the most famous BCS picture first introduced to de- 
scribe the phenomenon of superconductivity. There are 
many variants of the BCS theory available in the litera- 
ture. We will report here the predictions of the complete 
BCS scheme which accounts for non trivial renormaliza- 
tions of the physical quantities associated with quantum 
and thermal fluctuations. The many body solution pro- 
ceeds through a proper diagonalization of the Hamilto- 
nian (| 1 511 6[> by applying the Bogoliubov transformation 
to the Fermi field operators. This approach is non per- 
turbative and the ground state differs in a profound way 
from the uncorrelated wave function of the ideal Fermi 
gas, being characterized by peculiar correlations associ- 
ated with long range order. The ground state energy, 
expanded in terms of the small parameter kpa, takes the 
same form (fT5)l holding for the repulsive gas, pairing ef- 
fects being responsible only for higher order exponential 
corrections. Of course in this case the first correction, 
linear in a, gives a negative contribution to the ground 
state energy. BCS theory predicts a phase transition as- 
sociated with the occurrence of long range order. The 
corresponding critical temperature is given by the result 

T c = 0.28e 2lT /kFa T F (19) 

showing that the critical temperature becomes exponen- 
tially small as one decreases the value of £;_F|a|, mak- 
ing the observability of superfluid phenomena a diffi- 
cult task in dilute samples. Actually, in the experimen- 
tally relevant case case of harmonically trapped config- 
urations the predicted value for the critical temperature 
easily becomes even smaller than the oscillator temper- 
ature Huho/ks = Tp/(6N a ) 1 ^ 3 . For example, choosing 
kpa = —0.2, one finds T c ~ 10~ 4 TV a value significantly 
smaller than the oscillator temperature for realistic val- 
ues of N a . 

Thanks to the Fesbach resonances exhibited by several 
atomic species it is now possible to change the value of 
the scattering length in a highly controlled way by simply 
tuning the external magnetic field. For example, start- 
ing from a negative and small value of a it is possible 
to increase the size of the scattering length, reach the 
resonance where a diverges and explore the other side of 
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the resonance where the scattering length becomes pos- 
itive and eventually small. One would naively expect to 
reach in this way the regime of the dilute repulsive gas 
discussed above. This is not the case because in the pres- 
ence of the Feshbach resonance the positive value of the 
scattering length is associated with the emergence of a 
bound state and the formation of dimers. The size of 
the dimers is of the order of the scattering length a and 
their binding energy is ~ h 2 /ma 2 . These dimers have a 
bosonic nature, being composed of two fermions and if 
the gas is sufficiently dilute they consequently give rise to 
Bose-Einstein condensation, as experimentally proven in 
[20| . The size of the dimers cannot be however too small 
otherwise collisions between dimers give rise to transi- 
tions to deeper molecular states [19]. The behaviour of 
the dilute gas of dimers (kpa <C 1), hereafter called the 
BEC limit, is properly described by the theory of Bose- 
Einsten condensed gases available in both uniform and 
harmonically trapped configurations [3j, |g, |7fl . In partic- 
ular we can immediately evaluate the critical tempera- 
ture T c . In the uniform case this is given by the text- 
book relationship T c = (2nh 2 /M)(n a / 5 3 / 2 (l)) 2/3 where 
n a is the density of dimers (equal to the density of each 
spin species), M = 2m is the mass of the dimers and 
.93/2(1) = 2.612. In terms of the Fermi temperature of 
the uniform gas one can write 

T c = 0.22T F (20) 

showing that the critical temperature for the onset of 
superfluidity, associated with the Bose-Einstein conden- 
sation of dimers, takes place at temperatures of the or- 
der of the Fermi temperature, i.e. at temperatures much 
higher than the exponentially small value (fit?)) character- 
izing the BCS regime. For this reason one often speaks 
of high T c superfluidity. The inclusion of interactions be- 
tween dimers, fixed by the molecule-molecule scattering 
lenght dj\/ according to the relationship au = 0.6a [l9| . 
is also straightforward and is provided, at T = 0, by the 
Gross-Pitaevskii theory. 

The gas of dimers and the gas of repulsive fermions 
discussed above represent two different branches of the 
many-body problem, both corresponding to positive val- 
ues of the scattering length. The physical implementa- 
tion of the repulsive gas configuration can be achieved 
by switching on adiabatically the value of the scattering 
length starting from the value a — (see [3l|). Viceversa 
the gas of dimers is naturally implemented by crossing 
the Feshbach resonance starting from negative values of 
a (experimentally it is also realized by cooling down a 
gas at fixed and positive value of the scattering length) . 
The gas of dimers exploits the attractive nature of the 
force which is crucial in order to ensure the binding of 
the fermionic pairs. Furthermore, differently from the re- 
pulsive Fermi gas, at low temperature the gas of dimers 
exhibits Bose-Einstein condensation and superfluidity. 

A more difficult problem concerns the behavior of the 
many body system when the scattering length becomes 
larger than the interparticle distance. This corresponds 



to the challenging situation of a dilute (in the sense that 
the range of the force is much smaller than the interpar- 
ticle distance) but strongly interacting system. Will the 
system be stable or will it collapse as happens for bosons 
interacting with large scattering lengths? At present an 
exact solution of the many-body problem is not avail- 
able for kp\a\ > 1 and one has to make use of approx- 
imate schemes or numerical simulations. Our present 
understanding, based on both theoretical estimates and 
experimental results, is that the system remains stable 
even in the so called unitary limit, corresponding to 
kp\a\ 3> 1. Furthermore in this limit new interesting 
features are expected to take place. In fact the thermo- 
dynamic quanties should no longer depend on the actual 
value of the scattering length, the only relevant lenght 
scales of the problem being the inverse of the Fermi wave 
vector and the thermal wave length (2lj . An important 
consequence is that, at T — 0, the chemical potential can 
be parametrized in the simple way: 

M = (1 + 0)E F = (1 + /?)( 3 7r 2 ) 2 / 3 ^n 2 / 3 (21) 

2m 

where (3 is a universal dimensionless parameter (22I |23| 
and n = rtj+nj is the total density. The relationship (|2ip 
fixes the density dependence of the equation of state with 
non trivial consequence on the behaviour of the density 
profiles and of the collective frequencies of harmonically 
trapped superfluids (see the next sections). 

In general we have no reason to doubt that, at T = 

0, the system will be superfluid for all values of kpa, 

1. e. along the so called BEC-BCS crossover provided by 
the Feshbach resonance. This will be actually the basic 
point guiding our discussions in the following sections. 
The above results make also plausible to assume that 
for broad resonances, corresponding to kpR* <C 1, all 
the relevant properties of the system can be described 
in terms of the dimensionless combination kpa. This 
introduces a remarkable simplification in the theoretical 
description of this non trivial many-body problem. 

While the exact solution of the many-body problem 
along along the BCS-BEC crossover is not available, a 
useful approximation is provided by the BCS mean-field 
theory. This approach was first introduced to investi- 
gate the crossover by Eagles, Leggett and others [24[ 
with the main motivation to explore the properties of su- 
perconductivity and superfluidity beyond the weak cou- 
pling limit kp\a\ <C 1. The main merit of this approach 
is that it provides a comprehensive, although approxi- 
mate, description of the equation of state along the whole 
crossover, including the unitary limit 1 /kpa — > and the 
BEC regime of small and positive a. 

The idea of BCS mean field theory is based on the re- 
duction of the microscopic Hamiltonian to a mean field 
form where only the anomalous average < ^>^> > is taken 
into account, while the normal average < ty'ty > is ne- 
glected in order to ensure convergency to the resulting 
equations along the whole crossover. The anomalous av- 
erage characterizes the occurrence of long range order in 
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the off-diagonal 2-body density matrix and the non van- 
ishing value of the pairing function 

F(r,s)=<* T (r+|)tt ; (r-|)> (22) 

The diagonalization of the resulting Hamiltonian is ob- 
tained with the help of the Bogoliubov transformation 
which transforms particles into quasi-particles. This pro- 
cedure, applied to the regularized potential (I16p . yields 
the non trivial equation 



dk 



(2tt) 3 \h 2 k 2 



1 

2E~ k 



where 



E k = J A 2 
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(23) 



(24) 



is the energy of the elementary excitations fixed by the 
order parameter A and 



V 



h 2 k 2 

2m 



(25) 



is the energy of a free particle calculated with respect to 
the chemical potential. The order parameter is directly 
related to the pairing function (f2"2"| through the relation 



A 



dsV eff (s) < * T (r+-)* x (r- 



';) > (26) 



where V e j / is the pseudopotential ([16]) . For positive val- 
ues of the chemical potential the order parameter A co- 
incides with the energy gap of the particle excitation 
spectrum (|2"4")h For negative values of /it the gap is in- 
stead given by \/ A 2 + n 2 . The gap has been measured in 
trapped gases [25| along the BEC-BCS crossover through 
radio-frequency transitions. 

Equation (|23p provides an important relationship be- 
tween A and the chemical potential [i entering the single- 
particle energy rjk- A second relation is given by the 
normalization condition which is preserved by the Bo- 
goliubov transformation and takes the form: 



dk 



1 - 



E k 



(27) 



The two equations should be solved in a consistent way. 
It is remarkable that they admit a solution for any value 
of kpa. The resulting prediction for the equation of state 
n(n) given by the BCS-mean field theory is reported in 
Fig. 3 where it is compared with the ab initio calcula- 
tion of [26}, based on a Monte Carlo simulation, and the 
selfconsistent calculation of [27J based on a diagrammatic 
expansion. The figure reveals a general qualitative agree- 
ment between the two approaches. In particular in the 
BEC limit the mean field theory accounts for the exis- 
tence of bound molecules of energy h 2 /ma 2 . The BCS 
mean field theory however misses some important fea- 
tures which are worth being mentioned: 



i) in the BEC limit the interaction between molecules 
is known to be governed by the value am = 0.6a of the 
molecule- molecule scattering length . BCS theory in- 
stead provides the wrong value <2m = 2a 

ii) In the BCS limit the mean field BCS theory does not 
account for the leading corrections to the ground state en- 
ergy (see eq.([T8"|)). This is due to the fact that the theory 
ignores the Hartree terms < > in the calculation of 
the energy. 

iii) the value of the order parameter A in the BCS 
regime misses the Gorkov-Melik Barkhudarov correction 
which reduces the proportionality coefficient between A 
and the Fermi energy by a significant factor (~ 0.45). 

iv) At unitarity the BCS theory predicts the value (3 — 
—0.42 to be compared with the value (3 = —0.58 given 
by the numerical Monte Carlo simulation. 



IV. EQUILIBRIUM PROPERTIES OF A 
TRAPPED GAS 



Let us first point out that in the deep BEC regime, 
where the interacting Fermi gas behaves like a gas 
of weakly interacting dimers, systematic information is 
available from our advanced knowledge of the physics 
of dilute Bose gases in harmonic traps H). Although 
the BEC regime is not easily achieved in present ex- 
periments with ultracold Fermi gases, the correspond- 
ing predictions nevertheless provide a useful reference for 
important physical quantities. For example the critical 
temperature for BEC in a harmonic trap is given by the 

1 /3 

well known expression kBTsEC — 0.94/ia;^ o iV CT which, 
combined with the corresponding expression ([7]) for the 
Fermi energy, provides the useful relationship 



Tbec = 0.527> 



(28) 



Result (1281) reveals that, once expressed in units of the 
Fermi temperature, the critical temperature in the BEC 
regime is higher than the corresponding value (|20p hold- 
ing in uniform matter. In this sense superfluidity is 
"favoured" by the presence of the harmonic trap. The 
critical temperature in trapped Fermi gases can be esti- 
mated also in other regimes, using the results of uniform 
matter within local density approximation. First esti- 
mates of T c in the BCS regime were given in [28j . 

In the following we will show how the results of the 
previous section can be usefully employed to describe 
the equilibrium properties of a trapped gas along the 
BEC-BCS crossover at zero temperature. To this pur- 
pose we will make use of the local density approxima- 
tion (LDA). This approximation assumes that, locally, 
the system behaves like a uniform body so that the en- 
ergy of the trapped system can be written in the integral 
form 



E = J dr[e(n(r) + V ho (r)n(r)} 



(29) 
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given by the sum of the internal, also called release, en- 
ergy 



El- 



and of the oscillator energy 



dre(n(r)) 



E ho = J drV ho (r)n(r) 



(30) 



(31) 



provided by the trapping potential ([T]). In cas. (|29H3T|) 
n(r) = 7if(r) + nj_(r) is the total density profile deter- 
mined by the variational relation S(E — fi N)/8n(r) = 
which yields the most relevant LDA equation 



where 



Mo = K n ( r )) + Vho(r) 



(32) 



(33) 



is the chemical potential of uniform matter and (j,q is 
the chemical potential of the trapped gas, fixed by the 
normalization condition J dm(r) — N. Equation (|3"2")) 
provides an implicit equation for the density profile n(r) 
holding at equilibrium. 

The applicability of LDA in a Fermi gas is justified if 
the external potential is sufficiently smooth (which is in 
general guaranteed by the harmonic choice) and if the 
relevant energies are much larger than the typical os- 
cillator energies providing the quantization of the single 
particle levels (/io 3> tu^o)- In the absence of interac- 
tions the equation of state is given by the ideal Fermi gas 
expression 



ftfn) = (37T 2 ) 2 / 3 — n 2 / 3 
2m 



yielding the result 



n(*) = ^( 2 -^) 3/2 (»o-V ho (r) 3/2 



(34) 



(35) 



for the equilibrium density profile. For harmonic trap- 
ping this corresponds to the 3/2-th power of an inverted 
parabola (see eq.©). The chemical potential /x takes 
the value 



(36) 



and, as expected, coincides with the Fermi energy ([7| 
already introduced for the trapped ideal gas. 

Interactions modify the shape and the size of the den- 
sity profiles. The effects are accounted for by eq. ([32"|) 
once the equation of state fi(n) is known. A simple re- 
sult is obtained at unitarity where the equation of state 
(see eq. (l2"T10 has the same density dependence (|34f as for 
the ideal gas, apart from a dimensionelss renormalzation 
factor. By dividing (|3"2"|) by (1 + /3) one then finds that the 
results at unitarity are simply obtained from the ones of 



the ideal Fermi gas by a simple rescaling of the trapping 
frequencies and of the chemical potential. In particular 
the density profile at unitarity takes the same form 



-3/2 



as in the ideal gas, with the Thomas-Fermi radii given 
by the rescaled law 



(38) 



where Fq are the Thomas-Fermi radii of the ideal gas 
(see eq. |g2D). 

Another important case is the BEC limit where one 
treats the interaction between dimers using the mean 
field equation of state fi = guni '4 with the coupling con- 
stant gM = 4:irh 2 aM /2m fixed by the molecule-molecule 
scattering length. According to the well known results 
holding for weakly interacting bosons one finds 



Ri 



2 ~ 



aho 



(39) 



In Fig. 4 we show the experimental results [29| for the 
in situ density profiles in a harmonically trapped Fermi 
gas at extremely low temperature. The results actually 
correspond to the axial column density, i.e. to the in- 
tegral J dxdyn(r) which is the quantity measured in the 
experiment of [29j. Very good agreement between eperi- 
ment and theory is found at unitarity where the column 
integral of (|3"T|) is easily calculated and yields the result: 



_ 16 N _ z% 5/ 2 

n column( z ) — g (I r? 2 ' 



(40) 



with Z = R z given by ea. (|3"8")) The best fit to the ex- 
perimental curve yields the value [3 = —0.68 revealing 
explicitly the attractive role of interactions at unitarity. 
In Fig. 5 we show the experimental results for the axial 
radius of the trapped cloud along the crossover obtained 
in the same experiment. The results on the BEC side are 
reasonably consistent with the prediction (|39[) (continu- 
ous line) once one uses the value clm = 0.6a, provided by 
the exact solution of the dimer-dimer scattering problem 
(T^ |. For a systematic comparison between experimental 
and theoretical results of the density profiles see [3(| • 

In addition to the in situ density profiles a valuable 
source of information comes from the measurement of 
the release energy, i.e. the energy of the system measured 
after switching off the confining trap. The release energy 
is defined by the sum 



E r 



Ekin + E2- 



2 — body 



(41) 



of the kinetic and the 2-body interaction terms and, 
in the local density approximation, is simply given by 
eq. (f30"]) . i.e. by the first term of the total energy (|29j) . 
The integral (|30[) should be in general carried out nu- 
merically. There are however important cases where an- 
alytical results are available. For example, at unitarity, 
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one finds the useful result 
3 



E r 



:l 
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(1 + 13) 1/2 NE I 



(42) 



where Ep = huhoi^N) 1 ^ 3 is the Fermi energy (|36|) of the 
ideal trapped gas. 

From the measurement of the release energy T. Bour- 
del et al. [2(| were able to extract the value —0.64 for the 
parameter /3 in reasonably good agreement with the value 
extracted from the in situ measurement of the profile as 
well as with the ab initio predictions of theory. 

Another quantity that has been recently measured in 
ultracold Fermi gases is the momentum distribution [32| . 
To this purpose one switches off the scattering length, 
profiting of the existence of a Feshbach resonance, just 
before the expansion of the gas. The images of the 
density profile after expansion then provide direct ac- 
cess to the momentum distribution. The experimental 
and theoretical results are shown in Fig. 6 for different 
values of k F a. The theoretical curves have been ob- 
tained by applying the local density approximation [33| 
n(k) = (27r)~ 3 J drnfe(r) to the particle distribution func- 
tion 



z< a\cik >- 



1 - 



E K 



(43) 



calculated in uniform matter at the corresponding value 
of the density, within BCS mean field theory. Eq. (|43|) 
reduces to the step function 8(fc — fop) in the deep BCS 
limit, while it approaches the value 



(44) 



in the opposite BEC limit where, as expected, it coin- 
cides with the momentum distribution of a particle in 
the bound state of the dimer. The figure reveals a rather 
satisfactory agreement between theory and experiments. 



V. DYNAMICS AND SUPERFLUIDITY 

Superfluidity is one of the most challenging properties 
exhibited by ultracold Fermi gases, analog to the super- 
conductor behaviour taking place in charged Fermi sys- 
tems. It shows up in peculiar transport features. Among 
the most noticable manifestations one should recall the 
absence of viscosity, the hydrodynamic nature of macro- 
scopic dynamics even at zero temperature, the existence 
of quantized vortices, the occurrence of pairing effects 
and the vanishing of the spin polarizability. The last two 
features are typical of Fermi superfluids, while the first 
ones characterize also the superfiuid behavior of Bose sys- 
tems. The possibility of exploring these phenomena in 
ultracold gases provides a unique opportunity to comple- 
ment our present knowledge of superfluidity in neutral 
Fermi systems, previously limited to liquid 3 He. In this 



chapter we will discuss the hydrodynamic behaviour ex- 
hibited by superfluids and its implications on the dynam- 
ics of trapped Fermi gases. We will also discuss the dy- 
namic behaviour at higher energies where pair breaking 
effects become important. The implications of superflu- 
idity on the rotational properties will be discussed in the 
following section. 

The macroscopic behaviour of a neutral superfiuid is 
governed by the equations of irrotational hydrodynamics. 
At zero temperature they consist of coupled and closed 
equations for the density and the velocity field. In fact, 
due to the absence of the normal component, the su- 
perfiuid density coincides with the total density and the 
superfiuid current with the total current. The equations 
take the form 



9 , \ 

—n + V-(nv) = 
ot 

for the density (equation of continuity) and 



d 

TO — V 

dt 



V | -mv 2 + n(n) + V hl 



(45) 



(46) 



for the velocity field (Euler equation) where fi(n) is the 
chemical potential per particle, fixed by the equation of 
state of bulk matter. At equilibrium (v=0) the Euler 
equation provides the Thomas-Fermi condition (|32[) for 
the ground state profile. The irrotationality of the ve- 
locity field, implied by cq. (j4l)|) , is the consequence of the 
existence of the order parameter (|2"6"|) whose phase $ is 
related to the superfiuid velocity field by the relationship 



h 

2m 



(47) 



The hydrodynamic equations (|74H46|) differ from the cor- 
responding equations holding in the collisionless regime 
of a non superfiuid system because of the irrotationality 
constraint (|47|) . 

Despite the quantum origin underlying the superfiuid 
behaviour, the hydrodynamic equations of motion have a 
classical form and do not depend explicitly on the Planck 
constant. This peculiarity raises the question whether 
the hydrodynamic behaviour of a cold Fermi gas can be 
used to test the achivement of the superfiuid regime. As 
we will see, Fermi gases above the critical temperature 
can easily enter a collisional regime where the dynamic 
behaviour is governed by the same equations. In this 
respect it is important to stress that collsional hydrody- 
namics admits the possibility of rotational components in 
the velocity field which are strictly absent in the super- 
fluid. A distinction between classical and superfiuid hy- 
drodynamics is consequently possible only studying the 
rotational properties of the gas (see next section) . 

The applicability of the hydrodynamic equations is re- 
stricted to the study of macroscopic phenomena, charac- 
terized by low long wave length excitations. In particular 
the wave lengths should be larger than the so called heal- 
ing length. In the limit of BEC dimers the healing length 
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is proportional to y/h 2 /Mg^n where M — 2m and gun 
is the molecular mean field energy. In the opposite BCS 
limit the healing length is instead proportional to fij A 
where A is the pairing gap. At unitarity the healing 
length is of the order of the inter-particle distance. At 
the end of the section we will relate the healing length 
to the critical Landau's velocity and discuss its behav- 
ior along the BEC-BCS crossover. Let us finally remark 
that the hydrodynamic equations of superfluids have the 
same form both for Bose and Fermi systems, the effects of 
statistics entering only the form of the equation of state 
(j,(n). 

Important applications of the hydrodynamic equations 
concern the expansion of the gas after release of the trap 
and the collective oscillations. 

In most experiments with ultracold atomic gases im- 
ages are taken after expansion of the cloud. In the ab- 
sence of interactions the expansion of a Fermi gas is 
asymptotically isotropic even if the gas is initially con- 
fined by an anisotropic potential. This is the consequence 
of the isotropy of the momentum distribution n(p) which, 
for a non interacting gas, fixes the asymptotic beavior of 
the density n(r,t) of the expanded gas according to the 
law n(r, f) — > (m/i) 3 n(p) whith r = tp/m. Deviations 
from isotropy are consequently an important indicator of 
the role of interactions. In the experiment of [35[ et al 
the first clear evidence of anisotropic expansion of an in- 
teracting ultracold Fermi gas was reported, opening an 
important debate in the scientific community aimed to 
understand the nature of these novel many-body configu- 
rations. Hydrodynamic theory has been extensively used 
in the past years to analyze the expansion of Bose Ein- 
stein condensed gases. More recently it was proposed as a 
theory for the expansion of a Fermi superfluid [361 ] . The 
hydrodynamic solutions are obtained starting from the 
equilibrium configuration, corresponding to a Thomas- 
Fermi profile, and then solving eq. (|46|) by setting Vho = 
for t > 0. For an important class of configurations the 
spatial dependence can be analytically inferred. In fact, 
if the chemical potential has the power law fj, oc n 1 de- 
pendence on the density, the Thomas-Fermi equilibrium 
profiles have the analytic form no oc (/xq — Vho) 1 ^ 1 and 
one can easily prove that the scaling ansatz 

x y z 

n((x,y,z,t) = (b x b v b z )~ n Q (—, — , — ) (48) 

O x Oy Z 

provides the exact solution for the expansion with the 
scaling parameters bi obeying the simple time dependent 
equation: 



h(b x b y b z )^ 







(49) 



Equation (|4"9"j) generalizes the scaling equations previ- 
ously introduced in the case of an interacting Bose gas 
(7 = 1) [iH . From the solutions of ea.(|49p one can easily 
calculate the aspect ratio as a function of time. For an 



this is defined as the ratio between the radial and axial 
radii. In terms of the scaling parameters bi it can be 
written as 



R±(t) b ± (t) u z 



Z(t) b z (t) LU ± 



(50) 



For an ideal gas the aspect ratio tends to unity, while 
the hydrodynamic equations yield an asymptotic value 
7^ 1. Furthermore hydrodynamics predicts a peculiar 
inversion of shape during the expansion caused by the 
hydrodynamic forces which are larger in the direction 
of larger density gradients. As a consequence an initial 
cigar shaped configuration is brought into a disk profile 
at large times and viceversa. One can easily estimate 
the typical time at which the inversion of shape takes 
place. For a highly elongated trap (ui± ^> uj z ) the axial 
radius is practically unchanged for short times since the 
relevant expansion time along the z-th axis is fixed by 
l/u! z ^S> l/u>±. Conversely the radial size increses fast, 
and, for oj±t 1 one expects R±(t) ~ R±(0)u>±t. One 
then finds that the aspect ratio is equal to unity when 

LO z t r*> 1. 

In Fig. 7 we show the experimental images [35[ of the 
expanding cloud of an ultracold Fermi gas taken at dif- 
ferent times close to unitarity. The figure clearly shows 
the inversion of the shape from cigar to disc predicted 
by hydrodynamic theory. In Fig. 8 the predictions for 
the aspect ratio given by eqs. (|49H50p at unitarity, 
where 7 = 2/3, are shown together with the experi- 
mental results of [35| ■ The configuration shown in these 
figures corresponds to an initial aspect ratio equal to 
R±/Z = 0.035. The comparison strongly supports the 
hydrodynamic nature of the expansion of these ultracold 
Fermi gases. The experiment was repeated at higher 
temperatures and found to exhibit a similar hydrody- 
namic behaviour even at temperatures of the order of 
the Fermi temperature, where the system cannot be su- 
perfluid. One then concludes that in the normal phase 
the system exhibits a collisional regime. This is especially 
plausible close to unitarity where the scattering length is 
very large. 

The collective oscillations of a superfluid gas provide a 
further relevant source of information. These oscillations 
can be studied by considering the linearized form of the 
time dependent HD equations (|74I46|) . corresponding to 
small oscillations n = no + Sn ex.p(—iu>t) of the density 
with respect to the equilibrium profile no, where uj is 
the frequency of the oscillation. The linearized equations 
take the form 



uj 2 5n = V • 



noV(g^n) 



the velocity field being fixed by the equation 



<9v 



on 



(51) 



(52) 



axially symmetric trap (lj x 



Let us first consider the case of isotropic harmonic 
trapping (lo x = cu y = lu z = uJho)- A general class of 
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divergency free (also called surface) solutions is available 
in this case. They are characterized by the velocity field 
v oc X7(r e Ye m ), satisfying the condition V ■ v = and 
corresponding to the behavior (du/dn)Sn oc r Yg m for 
the density variation. Using the identity (<9/i/<9n)Vno = 
— VV/i , holding for the density profile at equilibrium, it 
is immediate to find that these solutions obey the equa- 
tion 



on on 



The resulting dispersion law 



Lo(£) = V~£iO h , 



(53) 



(54) 



is independent of the form of the equation of state, as 
generally expected for the surface modes driven by an 
external force. This result provides a model independent 
characterization of the hydrodynamic regime. The re- 
sult in fact differs from the prediction ui(£) = £uuho of 
the ideal gas model, revealing the importance of inter- 
actions accounted for by the hydrodynamic description. 
Only in the dipole case {£ = 1), corresponding to the 
rigid oscillation of the center of mass, interactions do not 
affect the frequency of these modes. In addition to the 
surface modes an important solution predicted by the hy- 
drodynamic equations in isotropic harmonic traps is the 
£ — 0, m = breathing radial mode whose solution can 
be found in analytic form if the equation of state is poly- 
tropic (jU oc n 1 ). In this case the velocity field has the 
radial form v oc r and the frequency is equal to 



u(m = 0) = \/'Sj + 2uj ho 



(55) 



For 7 = 1 one recovers the well known BEC result y/bcuho 
[37l ] . while at unitarity one finds 2wf lo . It is worth stress- 
ing that the unitary result to — 2u>ho is not limited to 
small amplitude oscillations and keeps its validity beyond 
the hydrodynamic approximation [381 ]. 

In the case of axi-symmetric trapping (uj x = u> y = 
luj_ =/= lu z ) the third component fvm of angular momentum 
is still a good quantum number and one also finds simple 
solutions of eci. (|53p . The dipole modes, corresponding 
to the center of mass oscillation, have frequencies oj±_ for 
in = ±1 and u z for m — 0. The oscillations where the 
velocity field is linear in the spatial coordinates exhibit a 
richer structure. The solutions with m — ±2 and m = ±1 
are surface excitations of the form v oc V(cc ± iy) 2 and 
v oc V(x ± iy)z, with frequency given, respectively, by 



w(m = ±2) = V2i 



and 



Lu(m — ±1) = J L} ± + 



(56) 



(57) 



independent of the equation of state. The dispersion |56|) 
of the radial quadruple mode has been recently tested 
experimentally along the crossover [39| . 



The m — solutions instead depend on the equation 
of state. For a polytropic dependence of the chemical 
potential state (// oc n 7 ) the problem can be solved in an 
analytic way. The solutions are characterized by a veloc- 
ity field of the form v oc W[a(x 2 + y 2 ) + bz 2 ], resulting 
from the coupling between the £ = 2 and £ = modes 
caused by the deformation of the trap. The correspond- 
ing frequencies are given by [40j 

W 2 (m = 0) = i[2( 7 + lVi + ( 7 + 2K 2 ± 
^4(7 + 1)M + (7 + 2)M + 4( 7 2 - 3 7 - 



Equation (|58|) reduces to the one derived in [37| in the in- 
teracting Bose case (7 = 1), while at unitarity (7 = 2/3) 
it coincides with the HD prediction app lied to the isoen- 
tropic oscillations of the ideal gas [U H2, HI, HH . For 
elongated traps (ui z <C uj±) the two solutions (|55|) reduce 



to to = ^2(7 + T)lo± and to = y/{3j + 2)/(7 + 1)lo z . 

Let us discuss in more detail the behavior of the com- 
pression modes at unitarity where, for elongated traps, 
one finds the prediction ^lO/Sui^ and y/ 12 /5u> z for the 
radial and axial oscillations, respectively. Experimentally 
both the two modes have been investigated in ultracold 
Fermi gases [45, 46]. In Fig. 9 we show the recent exper- 
imental results for the compressional radial mode taken 
from (46j. The agreement between theory and experi- 
ment at unitarity is remarkable confirming our under- 
standing of the dynamic behaviour in this highly corre- 
lated regime where the scattering length is much larger 
than the interparticle distance and the system exhibits 
a universal behavior. It is also worth noticing that the 
damping of the oscillations is smallest near unitarity. 

When we move from unitarity the collective oscilla- 
tions exhibits other interesting features. Theory predicts 
that in the BEC regime (7 = 1) the frequencies of both 
the axial and radial modes are higher than at unitarity 
(equal to 2uj± and i l /5/2u> z respectievly) . Furthermore 
the first corrections with respect to the BEC prediction 
can be calculated analytically, by accounting for the first 
correction to the BEC equation of state /i = gn pro- 
duced by quantum fluctuations. This is the so called 
Lee-Huang- Yang (LHY) correction first derived in the 
framework of Bogoliubov theory of a uniform Bose gas 
. The resulting shifts in the collective frequencies can 
be calculated analytically by solving the hydrodynamic 
equations (|53|) and treating the LHY term in a perturba- 
tive way. In the case of the most relevant m = radial 
breathing mode in a cigar- like configuration (lu z <C lu±_) 
one finds the following result for the relative frequency 
shift 11 



8lo _ 1O50F 
~ ~ 256 



\J a 3 M n(0) 



(59) 



where to = 2lo^_ is the unperturbed value. The shift is 
positive reflecting the repulsive nature of the interaction 
between molecules As a consequence, the dispersion law, 
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when one moves from the BEC regime towards unitar- 
ity, exhibits a typical non monotonic behavior, since it 
first increases, as a consequence of the LHY effect, and 
eventually decreases to reach a lower value at unitarity 
[49l ]. In general the collective frequencies can be calcu- 
lated numerically along the whole crossover by solving 
the hydrodynamic equations once the equation of state 
is known. Fig. 9 shows the predictions obtained using the 
equation of state of the MC simulation [5(| (see also [HTI ] 
and of BCS mean field theory [12]. The MC equation of 
state accounts for the LHY effect while the mean field 
BCS theory misses it, providing a monotonic behavior 
for the compressional frequencies as one moves from the 
BEC regime to unitarity. The equation of state is in 
both cases consistent with the correct value ylO/3W-L at 
unitarity. The accurate measurements of the radial com- 
pression mode shown in Fig. 9 confirm the prediction of 
the MC simulation, providing an important test of the 
equation of state and the first observation of the LHY 
effect. 

The behavior of the collective frequencies on the BCS 
side of the resonance exhibits different features. Theoret- 
ically one expects that when the system reaches the BCS 
regime the frequencies should be the same as at unitar- 
ity, the equation of state being governed by the same 2/3 
power law density dependence. However, things behave 
differently in experiments. The observed behavior is not 
a surprise and can be qualitatively understood by notic- 
ing that when one moves towards the BCS regime the 
critical temperature and the pairing gap become smaller 
and smaller and soon reaches values of the order of the 
trapping oscillator frequencies. Under these conditions, 
even assuming that the system be at zero temperature, 
one looses superfluidity and the system is eventually ex- 
pected to behave like a dilute collisionless gas whose col- 
lective frequencies, apart from minor mean field correc- 
tions, should approach the higher values 2uj± and 2o> z , 
respectively, for the radial and axial modes. Experimen- 
tally this transition is observed for the radial mode where 
the relevant trapping frequency is higher and occurs at 
about kp\a\ ~ 1. It is also associated with a strong in- 
crease of the damping of the collective oscillation. 

As already pointed out, hydrodynamic theory de- 
scribes correctly only the low frequency oscillations of 
macroscopic nature, corresponding to sound waves in a 
uniform body. When one considers higher excitation 
energies the dynamic response should also include the 
breaking of pairs into two fermionic excitations. The 
general picture of the excitations produced by a density 
probe can then be summarized as follows (for simplicity 
we consider a uniform gas): at low frequency the system 
exhibits a gapless phononic branch whose slope is fixed 
by the sound velocity c and hence by the compressibility 
of the gas according to the equation 



mc = n 



d/j, 
dn 



(60) 



tinuum of excitations starting from a given threshold fre- 
quency, above which one can break pairs. The value of 
the threshold frequency depends on the value of the total 
momentum carried by the perturbation. A first estimate 
is provided by the BCS mean field theory which predicts 
the following result for the threshold 



2A for [i > and q < 2^2/i 



hhJth = 2\/ (q 2 /8m — fi) 2 + A 2 elsewhere (61) 

The interplay between phonon and pair breaking excita- 
tions gives rise to different scenarios along the crossover. 
In the BCS regime the threshold occurs at low frequencies 
and the phonon branch very soon reaches the continuum 
of single particle excitations. The behavior is quite differ- 
ent in the opposite BEC regime where the gapless phonon 
branch extends up to high frequencies. At large momenta 
this branch actually looses its phononic character and ap- 
proaches the dispersion q 2 /4m of a free molecule. In the 
deep BEC limit the gapless branch coincides with the Bo- 
goliubov spectrum of a dilute gas of bosonic molecules. 
At unitarity the system is expected to exhibit an inter- 
mediate behavior, the discretized branch surviving up to 
momenta of the order of the Fermi momentum. A de- 
tailed calculation of the excitation spectrum [53j], based 
on a proper time-dependent generalization of the mean 
field BCS theory, is provided in Fig. 10. 

The results for the excitation spectrum provide a use- 
ful insight on the superfluid behavior of the gas in terms 
of the Landau's criterion according to which a system 
cannot give rise to energy dissipation if its velocity, with 
respect to a container at rest, is smaller than the Lan- 
dau's critical velocity defined by the equation 



ftw q 



(62) 



where Hio q is the energy of an excitation carrying mo- 
mentum q. According to this criterion the ideal Fermi 
gas is not superfluid because of the absence of a thresh- 
old for the single particle excitations, yielding v cr = 0. 
The interacting Fermi gas of Fig. 10 is instead superfluid 
in all regimes. By inserting result (fBTj) for the threshold 
frequency into eq. (|62p one can calculate the critical value 
of v due to pair breaking. The result is 



[V&+t?-n] 



1/2 



(63) 



In the deep BCS limit \a\kp — > (corresponding to 
A <C eq. (|63|) approaches the exponentially small 
value v cr = A /pp. On the BEC side the value (|6"3")) 
instead becomes larger and larger and the relevant ex- 
citations giving rise to Landau's instability are no longer 
single particle excitations but phonons, and the critical 
velocity coincides with the sound velocity. A simple es- 
timate of the critical velocity along the whole crossover 
is then given by the expression 



At high frequency one expects the emergence of a con- 



fer = Min (c, v^) 



(64) 
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Remarkably, one sees that v cr has a maximum near uni- 
tarity (see Fig. 11), further confirming the robustness of 
superfluidity at unitarity. The large value of v cr should 
show up in a visible reduction of dissipation in exper- 
iments where one moves an external impurity in the 
medium at tunable velocities. Experiments of this type 
have been already_performed in the case of Bosc-Einstein 
condensed gases |55j . 

An important physical quantity, directly related to the 
critical velocity, is the healing length defined as 



(65) 



Apart from a trivial numerical factor it coincides with 
the usual definition hj \JMgMn of the healing length on 
the BEC side and with the size of Cooper pairs in the op- 
posite BCS limit. The healing length provides the typical 
length scale above which the dynamic description of the 
system is safely described by the hydrodynamic picture. 
It is smallest near unitarity. 

The knowledge of the excitation spectrum and of the 
corresponding matrix elements of the density operator 
allows one to calculate the dynamic structure factor Q 



(66) 

where Kk and hui are the momentum and energy, re- 
spectively, transferred by the probe to the sample, 
dp is the fluctuation of the Fourier component pk — 



J^ exp— ik- Yj of the density operator, 



E n )/h are the usual Bohr frequencies and Q is the par- 
tition function. The definition of the dynamic structure 
factor is immediately generalized to other excitation op- 
erators like, for example, the spin density operator. In 
dilute gases the dynamic structure factor can be mea- 
sured with Bragg scattering experiments. 

The main features of the dynamic structure factor are 
best understood in uniform matter where the excitations 
are described in terms of their momentum. From the pre- 
vious discussion on the excitation spectrum one expects 
that, for sufficiently small momenta, the dynamic struc- 
ture factor be characterized by a discretized peak and by 
a continuum of single-particle excitations at higher en- 
ergy. At higher momentum transfer the behavior will de- 
pend crucially on the regime considered along the BEC- 
BCS crossover. In Fig. 12 we report the T = predictions 
for the dynamic structure factor S(k, uj) at relatively high 
wave vectors (fc = Akp) where the discretized branch is 
available only on the BEC side of the resonance. For such 
values of k the calculation of the dynamic response fac- 
tor can be usefully applied, within a LDA procedure, to 
estimate the response of the system in a trapped con- 
figuration. The results presented in these figures are 
based on a time-dependent generalization of the BCS 
mean field theory. This theoretical approach accounts 
for both phononic and single-particle parts of the excita- 



tion spectrum as well as for the corresponding hybridiza- 
tion phenomena. On the BEC side of the resonance one 
clearly sees a discretized peak correponding to the free 
molecule excitation energy h 2 k 2 /4m. It is remarkable to 
see that even at unitarity, where moleculae do not exist 
as independent excitations and the discretized peak has 
merged into the continuum of single particle excitations, 
the dynamic structure factor exhibits a pronounced peak 
at ~ h 2 k 2 /4m. On the BCS side of the resonance the 
molecular signatures are instead completely lost and the 
response is very similar to the one of an ideal Fermi gas. 
In the same figure we also show the magnetic structure 
factor which is obtained replacing the density operator 
Pi = Pi + Pi with the spin density operator pj — p± 
in eq. (|66[) . The magnetic operator does not excite the 
phonon mode and its strength is restricted to the contin- 
uum of single particle excitations. 

From the knowledge of the dynamic structure factor 
one can evaluate the static structure factor, given by its 
frequency integral 



h f°° 1 
S(k) = - S(k,u;)du, = -<J2< 



zk- (r 



> (67) 



S{k,u J ) = Q- 1 Y j e- pEmn \<^Pk\n>\ 2 8{hw-hu Jmn ) JJj 

m,n 



where, in deriving the last equality, we have used defi- 
nition (|6"6")) and the completeness relationship J2 n \ n >< 
n| = 1. The static structure factor is related to the two- 
body correlation function g(r) by the relationship 



S(k) = l + n dv[g{r) - l]e 



(68) 



Its behavior, evaluated at different points along the BEC- 
BCS crossover, is shown in Fig. 13 where we report both 
the calculation of the dynamic mean field approach [53| , 
which has been obtained by directly integrating the dy- 
namic structure factor, and the results of the ab initio 
Monte Carlo calculations of [111, obtained by Fourier 
transforming the pair correlation function g{r). The 
static structure factor decreases linearly at small k as 
a consequence of the phononic nature of the excita- 
tion spectrum, while it approaches the incoherent value 
S(k) = 1 for large wave vectors where only the terms 
i = j in the sum of eq. (|67p survive. It worth notic- 
ing that at intermediate values of k the static structure 
factor exhibits a maximum which, in the BEC regime 
kpa << 1, takes the value S(k) = 2 over an extended re- 
gion of values of k. The origin of this plateau is directlly 
related to the molecular nature of elementary excitations. 
In fact, if k is larger than the Fermi momentum, but 
still small compared to the inverse of the size a of the 
molecules, the probe mainly excites free molecules with 
energy h 2 k 2 /4m. Using the model independent f-sum 
rule result 



dujS(k. uj)lu = N 



h 2 k 2 

2m 



(69) 



and assuming, following Feynman, that a single 
molecular-like excitation with energy h 2 k 2 /4m exhausts 
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the integral, one finds the value S(k) = 2. Although this 
result holds only in the deep BEC limit, the figure shows 
that the enhancement of S(k) with respect to the inco- 
herent atomic value S(k) = 1 is clearly visible also at 
unitarity. 



VI. ROTATING FERMI GASES AND 
SUPERFLUIDITY 

Superfluidity shows up in spectactular rotational fea- 
tures. In fact a superfluid cannot rotate like a rigid body, 
due to the irrotationality constraint (|47|) imposed by the 
existence of the order parameter. At low angular ve- 
locity an important macroscopic consequence of irrota- 
tionality is the quenching of the moment of inertia. At 
higher angular velocities the superfluid can instead carry 
angular momentum via the formation of singular vor- 
tex lines. The circulation of these lines is quantized. In 
the presence of many vortex lines a regular vortex lat- 
tice is formed and the angular momentum acquired by 
the system takes the classical rigid body value. Both the 
quenching of the moment of inertia and the formation of 
vortex lines have been the object of fundamental investi- 
gation in the physics of quantum liquids and have been 
more recently explored also in dilute Bosc-Einstein con- 
densed gases. In this section we summarize some of the 
main features exhibited by dilute Fermi gases where im- 
portant experimental results are already available. We 
first discuss the macroscopic consequences of the irro- 
tationality constraint (moment of inertia and collective 
oscillations) and then some key features of quantized vor- 
tices. 

The moment of inertia O relative to the z — th axis is 
defined as the response of the system to a rotational field 
— SIL Z according to the relationship 



< l z >= no 



(70) 



where L z is the third component of angular momentum 
and the average is taken on the stationary configuration 
in the presence of the perturbation. In the limit of small 
angular velocity one can employ the formalism of linear 
response theory and write the moment of inertia in the 
form 



9 = 2Q 



< n\L z \m > | 2 



(71) 



Result ([72]) holds for both the Bose and Fermi ideal gas. 
It assumes lu x ^ co y , but admits a well defined limit when 
lj x — > Lu y . In the Fermi case, when the number of parti- 
cles is large, one can use the semiclasssical estimate for 
the radii yielding < x 2 >cx 1/uj 2 and < y 2 >cx 1/uj x . In 
this case eq. (|T2")) reduces to the rigid value of the moment 
of inertia: 



e, 



Nm < x + y > 



(73) 



(for a non interacting Bose-Einstein condensed gas, at 
T = 0, where the radii scale according to < x 2 >cx l/u> x 
and < y 2 >oc l/u> X) one instead finds that — > as 

LO x — ► UJy). 

Interactions can change the value of the moment of in- 
ertia of a Fermi gas in a profound way. The simplest way 
to calculate Q in the superfluid phase is to use the irrota- 
tional hydrodynamic equations developed in the previous 
section, by considering a trap rotating with angular ve- 
locity £1 and looking for the stationary solutions in the 
rotating frame. The resulting value for the angular mo- 
mentum < L z > then permits to evaluate the moment of 
inertia through definition (|70[) . A similar procedure was 
implemented experimentally to generate the rotation of 
a dilute Bose gas [58[ 

The equations of motion in the frame rotating with 
the trap are easily obtained by adding the term —£IL Z to 
the Hamiltonian. The hydrodynamic equations, in the 
rotating frame, take the form 



d 



V-(n(v-fi x r) =0 



(74) 



and 



m J^v + V Q mv2 + m(") + Vho - wiv • (n x r)^ = 

(75) 

where the rotating trap is now described by the time inde- 
pendent potential T4 Q . One sees that the rotation affects 
both the equation of continuity and the Euler equation. 
Here v = (fi,/2m)V$, where $ is the phase of the order 
parameter, is the superfluid velocity in the laboratory 
frame, expressed in terms of the coordinates of the ro- 
tating frame. In the presence of harmonic trapping an 
important class of stationary solutions is obtained mak- 
ing the ansatz 



v = aSJ(xy) 



(76) 



where \n > and E n are the eigenstates and eigenenergies 
of the unperturbed Hamiltonian and Q is the partition 
function. There is a simple case where the sum ([7T|) can 
be calculated explicitly. This is the ideal gas trapped 
by a deformed harmonic potential where the moment of 
inertia takes the analytic form (5?| 



e = 



mN 



■[<y 2 >-<x 2 >){ojI+cjI) 



2(oj 2 v <y 2 > -ul < x 2 >)] 



(72) 



The equilibrium density, derivable from eq. (|75[) by set- 
ting dv/dt — 0, has the same Thomas-Fermi form as in 
the absence of rotation, with renormalized values of the 
oscillator frequencies: 



- 2 



2 i 2 

uj x + a 

2 i 2 

uj y + a 



2aQ 
2afl 



(77) 



while the equation of continuity yields the relationship 



a = -SCI 



(78) 
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where 



< r 



x 2 > 



< y 2 + x 2 > 



(79) 



is the deformation of the atomic cloud in the x — y plane. 
For small angular velocities one finds S — e where 



(80) 



is the deformation of the trap. By evaluating the angular 
momentum L z = m J dr(r x v)n with the velocity field 
(j76[) one finally finds that the moment of inertia is given 
by the irrotational form 



6 = 5 2 O r 



11) 



which vanishes for an axi-symmetric configuration, point- 
ing out the crucial role played by superfluidity in the 
rotation of the gas. 

The measurement of the moment of inertia is not di- 
rectly accessible in dilute gases. However useful infor- 
mation on the rotational properties of the system can be 
obtained through the study of the collective oscillations. 
In this context a special role is played by the so called 
scissors mode, an oscillation of the system caused by the 
sudden rotation of a deformed trap. If the angle of rota- 
tion 0o of the trap is small compared to the deformation 
of the trap (|50|) the rotation of the potential |T]) produces 
the perturbation SV per t — {uj 2 —Lu 2 )<fioxy which naturally 
excites the quadrupole mode. The behavior of the result- 
ing oscillation depends in a crucial way on whether the 
system is normal or superfluid. In fact, while the restor- 
ing force is, in both cases, proportional to the square of 
the deformation parameter of the trapping potential, the 
mass parameter, being proportional to the moment of in- 
ertia, behaves quite differently in the two cases. In the 
superfluid the problem can be addressed by solving the 
hydrodynamic equations with the ansatz Sn cx xy and 
v oc Vxy. One easily finds that the oscillation around 
equilibrium is characterized by the frequency [59| 



to 



V2i 



(82) 



where ui\ = uj 2 + uj 2 . The result holds also for a tri- 
axially deformed trap where uj x ^ uj y ^ lu z . For a normal 
gas in the collisionless regime one finds two frequencies 
at lu± — \u> x ± u) y \. For a normal gas in the collisional 
regime one instead predicts an oscillation with the same 
frequency (|82p in addition to a low frequency mode of 
diffusive nature, caused by the viscosity of the fluid. The 
different behavior exhibited by the normal gas (both in 
the collisionless and collisional phases) reflects the rigid 
nature of the classical moment of inertia. The scissors 
mode, previously observed in a Bose-Einstein condensed 
gas [60], has been recently investigated also in ultracold 
Fermi gases [39| . At unitarity the experiment has con- 
firmed the correctness of the hydrodynamic frequency 
I , while deeply in the BCS regime the beating between 



the frequencies lu± — \uj x ±uj y \ has revealed the transition 
to the collisionless regime. If the gas is normal, but too 
deeply in the collisional hydrodynamic regime, the dif- 
fusive mode predicted by classical hydrodynamics would 
be located at too low frequencies to be observable in ex- 
periments and the scissors mode would look the same as 
in the superfluid phase. This is expected to be the case 
at unitarity just above the critical temperature, thereby 
making the distinction between the superfluid and the 
normal phase, based on the study of the scissors mode, 
a difficult task. 

More promising perspectives to distinguish between su- 
perfluid and collisional hydrodynamics are provided by 
the study of the collective oscillations excited on top of 
a rotating configuration. In fact in the presence of vor- 
ticity Vxj; the equations of collisonal hydrodynam- 
ics contain an addiditional term depending on the curl 
of the velocity field, which is absent in the irrotational 
equations of superfluid hydrodynamics. In the labora- 
tory frame the equation for the velocity field given by 
collisional hydrodynamics is given by 



d fl 
m —v + V ( -mv 2 + n{n) + V ho 



mvx(V x v)(83Q 



where we have omitted viscosity effects. The collec- 
tive oscilations correspond to linearized solutions with 
n = uq + 5n and v = Vo + Sv. In the presence of a steady 
rotation of the trap at angular velocity f2 the oscillation 
frequencies resulting from the equations of rotational hy- 
drodynamics differ from the ones of superfluid hydrody- 
namics since in the former case the steady velocity field 
Vo is given by the rigid body value f2 x r, while in the lat- 
ter case is given by the irrotational value — fWVxy. After 
generating the steady rotation of the gas by adiabatically 
ramping the angular velocity of a deformed trap, one can 
suddenly stop the rotation and the system starts oscil- 
lating performing rotations around the symmetry axis 
of the trap. In the superfluid this procedure will excite 
the scissors mode with frequency y/2uj_\_. In the case of 
rotational hydrodynamcis the behavior will instead be 
different, the scissors mode being coupled with the rigid 
rotation of the cloud. Under the condition £1 ^ e 2 ui± the 
resulting oscillation is characterized by the beating law 
(f)(t) = (fl/\/2Lj±)sin(V2LJ±t)cos(flt) [40[. 

The possibility of distinguishing between superfluid 
and rotational hydrodynamics is a unique opportunity 
provided by ultracold Fermi gases. In fact Bose-Einstein 
condensed gases, above T c , are usually extremely dilute 
and collisionless. Viceversa in the Fermi case the normal 
gas can be easily dominated by collisions even at tem- 
peratures smaller than the Fermi energy as proven by 
the behavior of the aspect ratio during the expansion of 
the unitary gas (see the previous Section). 

Let us now discuss the behavior of the superfluid Fermi 
gas at higher angular velocities where quantized vortices 
are formed. Recent experiments have confirmed their ex- 
istence along the BEC-BCS crossover (see Fig. 14). In 
these experiments vortices are produced by spinning the 
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atomic cloud with a laser beam and are observed after ex- 
pansion by ramping the value of the scattering length to 
positve values in order to increase their visibility. Quan- 
tized vortices were previously extensively investigated 
with Bose-Einstein condensed gases [6l| . Quantized vor- 
tices emerge as stable configurations if the angular ve- 
locity exceeds a critical value fixed by the energy cost 
needed for their production. 

A quantized vortex along the z axis is associated with 
the appearence of a phase in the order parameter 
given by the form expicjj where 4> is the angle around the 
z — axis. This yields the complex form 



A(r) = A(r±,z) exp(z<; 



(84) 



for the order parameter A where, for simplicity, we have 
assumed that the system exhibits axial symmetry and we 
have used cylindrical coordinates. The velocity field v = 
(H/2m)'V4> of the vortex configuration has a tangential 
form with modulus 



(85) 



2mr i 



which increases as one approaches the vortex line, in con- 
trast to the rigid body dependence v = Jl x r character- 
izing the rotation of a normal fluid. The circulation is 
quantized according to the rule 



v-dt 



rn 



(86) 



which is smaller by a factor 2 with respect to the case 
of a Bose superfluid with the same atomic mass. The 
value of the circulation is independent of the radius of 
the contour. This is a consequence of the fact that the 
vorticity is concentrated on the z-axis according to the 
law 



Vxv = 7T— <5 (2) (r ± )z 



(87) 



where z is the unit vector in the z-th direction. It deeply 
differs from the uniform vorticity V x v = 2f2 of the rigid 
body rotation. 

The angular momentum carried by the vortex is given 
by the expression 

drr x vn(r) = TV- 
holding if the vortex line coincides with the symmetry 
axis of the density profile. If the vortex is displaced to- 
wards the perifery of a trapped gas the angular momen- 
tum takes a smaller value. In this case the axial symme- 
try of the problem is lost and the order parameter cannot 
be written in the form (|8"4l . 

A first estimate of the energy of the vortex line is ob- 
tained employing macroscopic arguments based on hy- 
drodynamics and considering, for simplicity, a gas con- 
fined in cylinder of radial size R. The energy E v acquired 



by the vortex is mainly determined by the hydrodynamic 
kinetic energy {m/2) J drv 2 n which, employing the veloc- 
ity field (|85|) . yields the following estimate for the vortex 
energy: 



h R 

E v = N An— 

AmR 2 £ 



(89) 



where we have introduced the radius £ of the core of the 
vortex which fixes the distance below which the hydrody- 
namic expression for the kinetic energy no longer applies. 
This size is identified with the healing length whose value 
varies significantly along the BEC-BCS crossover, being 
particularly large in the BCS limit. Equation (|8"9"|) can be 
used to evaluate the critical angular velocity f2 c for the 
existence of an energetically stable vortex line. This value 
is obtained by imposing that the change in the energy 
E — Q C L Z acquired by the system in the frame rotating 
with angular velocity fl c be equal to E v . One finds 



fi B = 



h , r 

rlTl— 



2mR 2 £ 



(90) 



Applying this estimate to a harmonically trapped con- 
figuration with Rtf ~ R and neglecting the logaritmic 
term which provides only a correction of order of unity, 
we find fl c /ui± ~ hu>±/Ef lo where ui± is the radial fre- 
quency of the harmonic potential and Eho ~ rnu)\Wj, F is 
the harmonic oscillator energy of the trapped gas. The 
above estimate shows that in the Thomas-Fermi regime, 
where Eho 3> huj± , the critical frequency is much smaller 
than the radial trapping frequency, thereby suggesting 
that vortices should be easily produced in slowing ro- 
tating traps. This conclusion, however, does not take 
into account the fact that the nucleation of vortices is 
strongly inhibited at low angular velocities by the occur- 
rence of a barrier. For example in rotating Bose-Einstein 
condensates it has been experimentally shown that it is 
possible to increase the angular velocity of the trap up 
to values significantly higher than f2 c without generat- 
ing vortical states. Under these conditions the response 
of the superfluid is governed by the equations of irrota- 
tional hydrodynamics. 

A challenging problem concerns the visibility of the 
vortex lines. Due to the smallness of the healing length 
they cannot be observed in situ, but only after expansion. 
In particular the healing length is very small in the most 
interesting unitary regime. Another difficulty in reveal- 
ing the vortex lines is the reduced contrast in the den- 
sity with respect to the case of Bose-Einstein condensed 
gases. Actually, while the order parameter vanishes on 
the vortex line the density does not, unless one works in 
the deep BEC regime. Microscopic calculations of the 
vortex structure along the BEC-BCS crossover [63[ are 
in most cases based on the generalization of the BCS 
mean field theory to include non uniform configurations. 
In Fig. 15 we report the predictions obtained in [63| for 
the density profile and for the order parameter A as a 
function of the distance from the vortical line. The fig- 
ure clearly shows that the contrast in the density profile 
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becomes weaker and weaker as one approaches the BCS 
regime where, differently from the opposite BEC regime, 
the density is not affected by the presence of the vortex. 
Conversely the order parameter A always vanishes close 
to the vortex line. 

At higher angular velocities more vortices can be 
formed giving rise to a regular vortex lattice. In this 
limit the angular momentum acquired by the system 
takes the classical rigid body value and the rotation will 
look similar to the one of a rigid body, characterized by 
the law Vxv = 2ft. Using result (|57|) and averaging 
the vorticity over several vortex lines one finds the result 
Vxv = {h/2m)n v i, where n v is the number of vortices 
per unit area, so that the density of vortices is related to 
the angular velocity f2 by the useful relation 



Tin 



(91) 



which turns out to be a factor 2 larger than in the case of 
Bose superfluids with the same value of the atomic mass. 
Equation (f9"Tj) shows that the distance between vortices 
(proportional to 1/ y/n^) depends on the angular velocity 
but not on the density of the gas. In other words the 
vortices form a regular lattice even if the average density 
is not uniform as happens in the presence of harmonic 
trapping. This intriguing feature, already pointed out 
in the case of Bose-Einstein condensed gases, has been 
confirmed by the recent experiments on Fermi gases (see 
fig.14). 

The vortex lattice is responsible for an important bulge 
effect associated with the increase of the radial size of 
the cloud. In fact in the presence of a rigid rotation 
the effective potential felt by the atoms is given by 
Vho — (rn/2)S! 2 r^, corresponding to a modified equilib- 
rium density, whose Thomas-Fermi radii satisfy the new 
relationship 



n 2 



(92) 



showing, by the way, that at equilibrium the angular ve- 
locity cannot overcome the radial trapping frequency. In 
experiments where the vortex lattice is not formed in the 
presence of a stationary rotating trap this formula can be 
used to evaluate the effective value of the angular velocity 
by just measuring the aspect ratio. 

Important consequences of the vortex lines concern 
also the frequency of the collective oscillations. For ex- 
ample, using a sum rule approach [&4l ] it is possible to 
show that the splitting between the m = ±2 quadrupole 



frequencies is given by the formula 



w(m = +2) - u)(m = -2) = 2- 



(93) 



where t z = (L z ) /N is the angular momentum per particle 
carried by the vortical configuration. For a single vortex 
line £ z is equal to H/2, while for a vortex lattice i z is 
given by the rigid body value Sim < r\ >. In the latter 
case one recovers exactly the splitting 251 predicted by 
the equations of rotational hydrodynamics which, in the 
case of axi-symmetric configurations, yield the result [4(| 



u)(m = ±2) 



2w 2 ± 



n 2 ±n 



(94) 



for the frequencies of the two m = ±2 quadrupole modes. 
The equations of rotational hydrodynamics actually pro- 
vide the correct description of the collective oscillations 
of a gas containing a vortex lattice which, from a macro- 
scopic point of view, behaves like a classical gas rotating 
in a rigid way. 

The experimental production and measurement of 
quantized vortices has provided a definitive proof of su- 
perfluidity in these ultracold Fermi gases. Actually the 
esixtence of vortices as stable configurations can be re- 
garded as a proof of the absence of viscosity. In fact in 
the rotating frame the velocity field of a vortical configu- 
ration is not vanishing and any tiny presence of viscosity 
would bring the system into a rigid rotation, correspond- 
ing to a vanishing velocity in the rotating frame. 



A. Conclusions 

In this paper we have have summarized some features 
of the superfluid behavior exhibited by ultracold Fermi 
gases with special emphasis on the behavior of the col- 
lective oscillations. Other subjects related to the dynam- 
ics of these novel systems that would deserve careful in- 
vestigation concern Fermi gases in the presence of in- 
balance (iVf ^ JVj.) and Fermi gases of unequal masses 
(m T ^ mi). 
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FIG. 1: Release energy in units of the Fermi energy (0 as 
a function of the temperature calculated for an ideal Fermi 
(solid line), classical (dotted line) and Bose (dashed line) gas 
with the same number of atoms and the same mass. Here T® 
is the critical temperature for Bose- Einstein condensation for 
an ideal Bose gas. 
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FIG. 2: Evidence for quantum degeneracy effects in trapped Fermi gases. The average energy per particle, extracted from 
absorption images, is shown for two spin mixtures. In the quantum degenerate regime the data agree well with the ideal Fermi 
gas prediction (solid line). From [{J. 




FIG. 3: Equation of state along the BEC-BCS crossover as a funtion of the dimensioneless parameter 1/kpa. The results of 
the fixed node diffused Monte Carlo calculations (FN-DMC) of [2(| and of the diagrammatic expansion of [27| are compared 
with the prediction of the BCS mean field theory. 
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FIG. 4: Experimental results for the column density profiles along the BEC-BCS crossover for a gas of 6 Li atoms. The 
continuous curve at 850G, corresponding to unitarity, is the best fit based on eq. (|40l From [29|). 
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FIG. 5: Experimental results for the axial radius along the BEC-BCS crossover for a gas of Li atoms. The data are plotted 
after normalization to the non interacting Fermi gas. The full line corresponds to the prediction of Eq. (139[) with ajwr = 0.6a. 
From [H). 




FIG. 6: Theoretical (left) vs experimental (right) values of the momentum distribution in a gas of 40 K atoms along the 
BEC-BCS crossover. From 
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FIG. 7: Images of the expanding cloud of an ultracold Fermi gase of 6 Li atoms at different times. The figure clearly shows the 
typical inversion from the cigar to the disc shape revealing the hydrodynamic nature of the expansion. From [35| . 
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FIG. 8: Aspect ratio as a function of time during the expansion of an ultracold Fermi gas at unitarity (red points: experiment; 
red line: hydrodynamic theory). For comparison the figure also shows the results in the absence of interactions (blue points: 
experiment; blue line ballistic expansion). From [351 ]. 
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FIG. 11: Landau's critical velocity (in units of Fermi velocity) 
calculated along the crossover. The figure clearly shows that 
the critical velocity is largest close to unitarity. From [Hit ]. 
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FIG. 12: Dynamic structure factor of the superfiuid Fermi 
gas along the BEC-BCS crossover at k — 4fc_F- From [53| . 




FIG. 13: Static structure factor of the superfiuid Fermi gas 
along the BEC-BCS crossover. From 
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FIG. 14: Experimental observation of quantized vortices in a 
superfiuid Fermi gas along the BEC-BCS crossover [62l ]. 



FIG. 15: Radial profiles of the density and of the order pa- 
rameter of a vortex line along the BEC-BCS crossover. From 
R. Sensarma et al. [63l ]. 



